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This report describes two projects. Firstly, a Reynolds stress closure for near- wall 
turbulence is described. It was motivated by the simpler k — e — v 2 model described 
in last year’s annual research brief. 

DNS of three-dimensional channel flow show a curious decrease of the turbulent 
kinetic energy. The second topic of this report is a model which reproduces this 
effect. That model is described and used to discuss the relevance of the three- 
dimensional channel flow simulation to swept wing boundary layers. 

1* Motivation and objectives 

The region of a flow very near to a surface exerts a disproportionate control on 
transport of momentum and heat between surface and fluid. The turbulence in this 
region is strongly inhomogeneous and highly anisotropic. The anisotropy is particu- 
larly important: the normal component of turbulent intensity (v 2 ) is responsible for 
transport to and from the surface; this component is suppressed by the proximate 
surface. The behavior of v 2 was the focus of attention in our previous work on 
near- wall modeling (Durbin 1991). An elliptic relaxation equation was formulated 
to describe the suppression of v 2 and to allow appropriate boundary conditions to 
be satisfied. Last year’s report described an application of the k — e — v 2 model to 
boundary layers and heat transfer. 

During the period covered by the present report the near-wall formulation was 
extended to a full Reynolds stress closure. This formulation is tensorally invariant: 
v 2 is not singled out as the normal component; the full Reynolds stress model is 
applicable to general geometries. The application of this model to separated flows 
is described in the article by S. Ko in the present volume. The model and an 
assessment of its predictive potential are described in section 2.1. 

Moin et al (1990) performed a DNS of three-dimensional channel flow in order 
to provide detailed data for use in modeling. This flow was meant as an idealization 
of the swept-wing boundary layer. We have attempted to model three-dimensional 
channel flow and to assess the relevance of this flow to the infinite swept wing 
boundary layer. The most peculiar feature observed in the DNS was a decrease 
of the turbulent intensity when the cross-stream pressure gradient was applied. 
One would expect the cross-stream flow to add a new contribution to turbulence 
production and to increase the turbulence intensity; indeed, this is what exist- 
ing turbulence models predict. To account for the reduction in intensity, we have 
assumed that the streamwise vorticity produced by the cross-flow enhances the cas- 
cade of energy to small scales. The mechanism might involve breakup of large, 
streamwise eddies. The present model consists simply of a parameterization of 
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three-dimensionality and an ad hoc modification to the c equation. With only this 
modification, the Reynolds stress model reproduces many of the features observed 
in the DNS. Three-dimensional shear flow is discussed in section 2.2. 


2. Accomplishments 

2.1 The Reynolds stress model 

A detailed description of the Reynolds stress model is provided in Durbin (1992a); 
the following is a brief summary. 

The exact Reynolds stress transport equation can be written 


DtuTu] = pij 4- Pa - e - dkUkUiUj + i/V 2 u,u, 

where D t (») is the convective derivative following the mean flow; 

Pij = —UiuidkUj - uJukdkUi 

is the rate of turbulence production by mean velocity gradients, and 


Pij 


1 l 5— UiUj 

= — uidjp Ujdip - Cij -\ —e 


(1) 

( 2 ) 
( 3 ) 


is the redistribution tensor. (The trace of (3) is not equal to zero, but it is common 
practice to ignore the trace of p , j on the grounds either that it is negligible or that 
it is a ‘pressure diffusion’ term that can be absorbed in the self-transport model.) 
The Rotta model for anisotropic dissipation, -Uiu]e/k, has been added to (1) and 
subtracted from (3): among other virtues, this has the effect of making pij vanish 
at rigid, no-slip boundaries. By definition, e = Su/2 is the rate of dissipation of 
turbulent kinetic energy ( k ). 

In the notation of equation (3), all unclosed terms have been incorporated into 
except for the transport terms. We follow the usual practice of modeling turbulent 
self- transport by gradient diffusion: 

—dituuju] = di (^^drnWh'j ■ ( 4 ) 

This amounts to regarding u,Uj as a ‘substance’ being transported by the turbulent 
velocity u t . The diffusional model is as much a representation of the smoothing 
effect of ensemble averaging as of convective transport per sei this is why the model 
is parabolic rather than hyperbolic. For the eddy viscosity, 

VTij = C^vTu] T (5) 


will be used. For the time scale T, we adopt 

T = max (—■, Cr(j) 7 )- 


( 6 ) 


Reynolds stress closure 


187 


This becomes k/e far from boundaries. Near a surface where k 0, it becomes the 
Kolmogorov time-scale Cx(i//e) 1 / 2 , which is a suitable lower bound on T. 

A second order closure for pij is simply a proposed relationship between this 
unknown and the dependent variable of equation (1), UiUj. In all closures to 
date, this relationship has consisted of algebraic formulae. However, those formulae 
are based on quasi- homogeneous assumptions — most notably in the rapid pressure 
strain term — which are incorrect in the strongly inhomogeneous near wall region. 
I have proposed in previous reports that elliptic effects within the flow that are 
caused by the proximity of a boundary might be included by formulating an elliptic 
relaxation model for pij. Because the exact elliptic effects are inherently non-local, 
they cannot appear explicitly in a single point model; they are represented quite 
indirectly by the present model. 

The elliptic relaxation model can be put into coordinate independent form and 
written 

pij = kfij (7) 

L 2 V 2 fij - fij = -n ij. (8) 

The length scale L is formulated by analogy to (6): 


L = Cl max 



( 9 ) 


Boundary conditions influence the solution in the interior of the flow through the 
homogeneous solutions to (8) (Durbin 1992a). Far from the surface, these solutions 
decay, and the f t j relax to quasi-homogeneity, as represented by a balance of the 
second term on the left side of (8) with its right side. 

Any quasi-homogeneous model can be used for II, j. Here we adopt the simple 
model recommended by Launder (1989). This consists of a sum of Rotta’s return 
to-isotropy and the isotropization-of-production models: 


Uii ~ kT* ( UiUj 3 kSiJ ) k ( Pii 3 PS,j ) (10) 

where 2 P — P t i, and Pij is given by (2). C\ is the usual Rotta constant, and 
C 2 =3/5 is found by the method of perturbation of isotropy. 

2.1.1 Assessment of the model 

The viability of the present Reynolds stress closure was assessed by comput- 
ing various boundary layer flows. An initial, detailed evaluation of the approach 
consisted of comparing the model solution to DNS data on channel flow; this com- 
parison is provided in Figure 1. The near wall behavior of the components of the 
Reynolds stress tensor is captured very well by the model. 

Durbin (1992a) contains further comparisons to data in zero pressure gradient, 
adverse pressure gradient, and curved wall boundary layers. I will not repeat all 
of those comparisons here, but simply provide a few representative figures. I will 
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FIGURE 1. Comparison of model (lines) to DNS (s ymb ols) profiles of Reynolds 
stresses in channel flow at R r = 395. +, n 2 ; x, Ar; A, u> 2 ; o, t> 2 ; ■, — uv. 



FIGURE 2. Friction coefficient (•) and displacement thickness (x) versus down- 
stream distance for the Samuel and Joubert experiment. Curves are solution of 
model, symbols experimental data. 
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FIGURE 3. Mean velocity profiles at stations 9 (x = 1.87m, •) and 12 (x = 
2.55 m, A) of the Samuel and Joubert adverse pressure gradient boundary layer. 
Lines are model, symbols experiment. 



FIGURE 4. Reynolds shear stress profiles at stations Tl ( x = 0.19m, x, ), 

T4 (x = 1.53 m, A, ) and T6 (x = 2.54 m, + , ) of Samuel and Joubert 

(1974). Curves are model, symbols experimental data. 
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skip the zero pressure gradient case. It provides a necessary prerequisite for any 
near wall turbulence model. The present model gives an excellent solution for the 
dependence of skin friction on Reynolds number and quite satisfactory profiles of 
mean velocity and turbulent Reynolds stresses. 

The Samuel and Joubert (1974) experiment on a boundary layer progressing into 
a region of increasingly adverse pressure gradient is a good test of the response of a 
model to pressure gradients. Solutions for skin friction and displacement thickness 
versus downstream distance are shown in Figure 2 along with experimental data. 
Profiles of mean velocity and Reynolds shear stress are contained in figures 3 and 4. 
Agreement between the model and data is quite reasonable, although discrepancies 
clearly can be seen. It is encouraging that the model is able to describe the departure 
from equilibrium observed in the experiment; in particular, Figure 4 shows how the 
pressure gradient causes a change of the near wall region from a constant stress 
layer to a layer of increasing stress. The log-layer theory invokes a constant stress 
layer. 

A nice illustration of the capability of the Reynolds stress model is provided by 
the boundary layer on a convex curved surface. A tensorially consistent model is 
formally independent of coordinate system. This means that the model equations 
can be projected onto curvilinear coordinates by the methods of differential geom- 
etry. Upon doing so, metric terms arise. The stabilizing effect of convex curvature 
can be attributed to those metric terms. Thus, much of the important physics are 
retained by the coordinate invariant Reynolds stress model. By contrast, nonin- 
variant models (e. g., mixing length) can not account for curvature effects, nor can 
quasi-isotropic, scalar models (e. g., k — e). 

Figure 5 shows the downstream evolution of the skin friction coefficient in a 
boundary layer that starts on a flat plate, then flows around a 90° convex, circular 
arc, and then continues along another flat section. This boundary layer was studied 
experimentally by Simon et al. (1982). In the experiment, the wind tunnel wall 
opposite to the test surface was contoured to maintain a constant surface pressure 
beneath the boundary layer. A pressure gradient normal to the curved wall exists 
to balance the centrifugal acceleration. 

The dashed line in Figure 5 shows the downstream evolution of C/ on a flat 
plate. Convex curvature begins at x = 0; the skin friction drops abruptly. In 
Figure 5, Cj has been normalized by its value at the start of curvature. This 
normalization enables one to invoke Reynolds analogy and include data on heat 
transfer coefficients in the figure. The heat transfer data are more comprehensive 
than the skin friction data and they were measured directly; the skin friction had to 
be inferred by assuming that the mean velocity profiles conformed to the universal, 
constant stress log-law. 

In fact, in the curved section the Reynolds shear stress profile does not show a 
constant stress layer near the surface. Figure 6 shows how the centrifugal stabi- 
lization suppresses -uv in the outer region of the boundary layer. Indeed, -uv 
becomes slightly negative in the data denoted by x when y/S 99 > 0.4. After the 
curved section, -uv initially recovers near the surface, and the region of increased 
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FIGURE 5. Solid line is computed skin friction, normalized on its value at the 
start of curvature; dashed line is normalized skin friction computed for a flat plate 
boundary layer. Experimental data on skin friction (x) surface heat flux (•) from 
Simon et al. (1982). 
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Figure 6. Model solutions and experimental data on Reynolds shear stress. The 

profiles are at x = -0.062m (#, ), x = 0.162m (x, ) and x - 1.124m (A, 

). These stations are upstream of the bend, 20.6° around the bend and in 

the downstream recovery region. 
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stress then propagates across the boundary layer. The model shows some discrep- 
ancies with the data. These are most significant in the recovery region downstream 
of the bend. The model recovers to the flat plate skin friction too rapidly (Figure 


5). 

2.2 Effects of three-dimensionality 

The work described in this section is presented at greater length in Durbin 

(1992b). , . , . 

The most clear-cut feature of a three-dimensional boundary layer is the skewing 

of the direction of the mean velocity vector with height. This may alternatively 
be described as a presence of mean vorticity in the streamwise direction: in the 
boundary layer approximation, the vorticity vector is SI = (5 y W, 0, -3 v f/); the 
projection of the vorticity onto the mean velocity is 


n u = u 2 d v 


w 

u 


= U 2 dy tan(/?) 


( 11 ) 


where (3 is the direction of the mean velocity relative to the x-axis. In a three- 
dimensional boundary layer, the quantities in (11) are non-zero; hence, the vorticity 
and velocity are not orthogonal, and the direction of the mean flow skews with 

An additional feature of a three-dimensional turbulent shear layer is that the 
Reynolds stress is not collinear with the mean shear: 


d y W 

d y U 



uv 


dyUvW — dyWuV ^ 0. (12) 

This is the most obvious feature of the turbulence, and has been demonstrated in 

numerous experiments. 

Other effects of three-dimensionality on the turbulence are less clear cut. Many 
experiments on three-dimensional boundary layers use strong adverse pressure gra- 
dients to turn the flow. The role of three dimensionality in this type of experiment 
is clouded by the presence of the adverse pressure gradient. In order to isolate 
three- dimensional effects, and to provide comprehensive, accurate data, Mom et 
al. (1990) performed a DNS of pressure driven three-dimensional channel flow. This 
is a horizontally homogeneous flow produced by subjecting a fully developed two- 
dimensional channel flow to an accelerating span wise pressure gradient. 

Perhaps the most curious feature of this flow is that the turbulent intensity ini- 
tially decreases. The cross-stream acceleration adds a cross-stream component of 
mean shear, increasing the magnitude of shear, so one would anticipate an in- 
crease in the turbulent intensity. Indeed, this flow will ultimately evolve into a 
two-dimensional channel flow, at an angle to the original flow, with a higher mean 
velocity and higher turbulent intensity. The initial evolution— which is all that 
could be numerically simulated— is counter intuitive. 
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The objective of the present study is to model the behavior of three-dimensional 
channel flow and then to assess its relevance to the swept wing boundary layer. I 
considered the possibility that the decrease of the turbulent intensity was due to a 
suppression of energy redistribution into v 2 y as proposed by Moin et ai This was 
unable to explain the effect. The most direct mechanism, and that which I modeled, 
is to suppose that the three-dimensionality increases the rate of the energy cascade 
to small scales. 

A parameter is needed to characterize three-dimensionality. This ought to be 
related to the presence of streamwise vorticity (11) and to the misalignment of 
stress and shear (12). Consider projection of the production tensor on the direction 
of the mean vorticity: in a two-dimensional shear layer, P,jf2*/|fl| is zero, but not 
in a three-dimensional layer. This suggests that the invariant Q, /|S 7 | 2 might 
characterize three-dimensionality. Unfortunately, this quantity is identically zero 
in parallel shear flow. A non-vanishing invariant is P 2 D = where 

Pfj = PikPkj\ this will be used to characterize the effect of three- dimensionality on 
e. In parallel shear flow, 

= (dyUvw - dyWuv f 
so Pso is related to (12), as required. 

The e-equation both closes the turbulent kinetic energy ( k ) budget and provides 
a time-scale for use in modeling: T = k/e\ or in viscous regions, T = (i^/e) 1 / 2 . The 
exact equation for e is of little value to modeling; the usual model equation is based 
on an ad hoc notion that the production and dissipation of e can be represented by 
a function of the production and dissipation of energy. Thus, the evolution equation 
for homogeneous turbulence is assumed to be of the form 

i = fHP/e) ( 13 ) 

where P = P„/2 is the rate of energy production and F(#) is some function. ( 13 ) 
is simply a dimensionally consistent form for expressing the assumed dependence 
on production and dissipation of energy. The standard e-equation is obtained by 
letting F(#) be linear: 

F(P/e) = C Ci P/e — C €2 ( 14 ) 

where the C € ’s are empirical constants. Here we simply extend the argument list 
to include P^: 

i = ^F(P/e,P 2 D /e 2 ). ( 15 ) 

Again, we adopt a linear model 

• P ~ Cc 2 € + C C 3 P 2 D /e 

e = . ( 16 ) 


The value of C t 3 = 4 was chosen. 
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x (m) 

FIGURE 9. Skin friction coefficient and direction of surface stress, in radians, 
for infinite swept wing boundary layer. Data of Elsenaar and Boelsma (1974): x, 

+; data of Pontikos and Bradshaw(1985): A, y. Model: , C t3 — 4.0; , 

C ea = 0. 



FIGURE 10. Kinetic energy (upper curves) and t> 2 for infinite swept wing; data 

from Elsenaar and Boelsma (1974). (•, x), x = 0.3m; (A, y), x = 0.6m; 

■■■, x = 0.6m, = 0. Normalization is by friction velocity at x = 0 and by 99% 

boundary layer thickness. 
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The model described in §2.1 was solved in conjunction with (16). The evolution 
of the turbulent kinetic energy with time is displayed in Figure 7 along with DNS 
data. It is seen that the model produces the initial drop of kinetic energy. Figure 
8 shows the behavior of individual components of the Reynolds stress tensor. They 
evolve with time in much the same manner as observed in the Moin et at. DNS. 

The model was solved for the infinite swept wing flow studied by Elsenaar and 
Boelsma (1974) and by Bradshaw and Pontikos (1985). This flow was produced by 
subjecting a boundary layer to a pressure gradient directed at 35° the initial flow 
direction. The boundary layer originated either at a swept leading edge (Elsenaar 
and Boelsma) or at a swept suction slot (Bradshaw and Pontikos) so that the 
whole flow was invariant with respect to translation along lines of constant pressure 
(ignoring end walls). 

The upper curves in Figure 9 show the skin friction coefficient versus downstream 
distance and the lower curves show the angle of the surface stress to its upstream 
direction. The solid lines are a solution with C e , = 4. The dashed lines are a 
solution with C t3 = 0, so that the 3D effect found in the channel flow simula- 
tions can be assessed. The two sets of experimental data show how poorly repro- 
ducible this experiment is. Given this ambiguity and an ambiguity in the pressure 
gradient prescribed in the present computation, it cannot be concluded that the 
three-dimensional effect shown by Figure 9 is significant. 

Figure 10 shows profiles of k and v 2 at two downstream locations. The first is 
near the beginning of the pressure gradient and the second is well into the three- 
dimensional region. The sets of curves at x = 0.6 m show how the model predicts 
that the 3D effect occurs near to the surface, where it causes a reduction of k. The 
experimental data suggest that this effect might be present in the flow. 

My conclusion is that the three-dimensional effects uncovered in the channel flow 
DNS have only a minimal influence on the swept wing boundary layer. This is 
probably because the cross-stream pressure gradient in the DNS was quite large 
(10 times that of the initial equilibrium flow). In the swept wing experiments the 
perturbation to the boundary layer was not so abrupt. 

3. Future plans 

The Reynolds stress model is being solved in increasingly complex flows. The 
article by S. Ko in this volume describes computations of separated boundary layers. 
This work on complex flows will continue. A code for computing flow in curvilinear 
geometries is being developed (with N. Mansour). This will be used to compute 
turbulent flow on aerofoils. 

The possibility of formulating an eddy viscosity transport model for use in com- 
plex flows is being explored. Although this is far less satisfactory than Reynolds 
stress modeling, simple models are preferred by those concerned with complex aero- 
dynamic flows. 
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